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Based on general symmetry considerations we investigate how the dependence of the tunneling 
anisotropic magnetoresistance (TAMR) on the magnetization direction is determined by the specific 
form of the spin-orbit coupling field. By extending a phenomenological model, previously proposed 
for explaining the main trends of the TAMR in (001) ferromagnet /semiconductor/normal-metal 
magnetic tunnel junctions (MTJs) [J. Moser et al, Phys. Rev. Lett. 99, 056601 (2007)], we 
provide a unified qualitative description of the TAMR in MTJs with different growth directions. In 
particular, we predict the forms of the angular dependence of the TAMR in (001), (110), and (111) 
MTJs with structure inversion asymmetry and/or bulk inversion asymmetry. The effects of in-plane 
uniaxial strain on the TAMR are also investigated. 
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I. INTRODUCTION 

The tunneling anisotropic magnetoresistance (TAMR) 
effect refers to the dependence of the magnetoresis- 
tance of magnetic tunnel junctions (MTJs) on the ab- 
solute orientation(s) of the magnetization(s) in the fer- 
romagnetic lead(s) with respect to the crystallographic 
axes.— >^>^>^ Unlike the conventional tunneling magnetore- 
sistance (TMR) effect, the TAMR is not only present 
in MTJs in which both electrodes are ferromagnetic but 
may also appear in tunneling structures with a single 
magnetic electrodeiii^ Because of this remarkable prop- 
erty, if the major challenge of increasing the size of the 
effect at room temperature is solved, the TAMR could 
be an attractive candidate for applications in the design 
of new spin- valve based devices whose components could 
operate with a single magnetic lead. In what follows we 
focus our discussion on the case of MTJs in which only 
one of the electrodes is ferromagnetic. 

The TAMR has been experimentally and theoreti- 
cally investigated in a variety of systems under dif- 
ferent configurations.i^'^-i'^'^i'^^^'ii^i^^i^^ii^ 
This diversity has made it difhcult to build a unified the- 
ory of the TAMR. In fact, although there exists a general 
consensus in identifying the spin-orbit coupling (SOC) as 
the mechanism responsible for the TAMR, it has been 
recognized that the way the SOC influences the TAMR 
may depend on the considered system and configuration. 

Two different configurations, the in-plane and out-of- 
plane configurations, have been considered for investigat- 
ing the TAMR (for an extensive discussion see Ref. [isl 
and references therein). The in-plane TAMR refers to 
the changes in the tunneling magnetoresistance when the 
magnetization direction, defined with respect to a fixed 
reference axis [x] , is rotated in the plane of the ferromag- 
netic layer. The in-plane TAMR ratio is defined asi^ 
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where R{9, 0) denotes the tunneling magnetoresistance 
for the magnetization oriented along the direction defined 



FIG. 1: Schematics of a MTJ composed of a normal-metal 
lead (bottom layer) , a semiconductor or insulator barrier (cen- 
tral layer), and a ferromagnetic electrode (top layer). The 
vector m indicates the magnetization orientation, while [x\ 
denotes a reference crystallographic axis. 



by the unit vector m = (sin 6' cos 0, sin 6' sine/), cos 6*) (see 
Fig. [J). 

In the out-of-plane configuration, the TAMR measures 
the changes in the tunneling magnetoresistance when the 
magnetization is rotated within the plane defined by the 
reference axis [x] and the direction normal to the ferro- 
magnetic layer. The out-of-plane TAMR is given by*i^ 

TAMR[,j {6) _ . (2) 

An important property of the TAMR is the form of its 
angular dependence. It has been experimentally shown 
that both the in-plane and out-of-plane TAMR exhibit 
rather regular, and relatively simple angular dependence 
with a well defined symmetry, in spite of the highly com- 
plicated band structure of the considered systems j^i^Siii^ 
This suggests that although the size of the TAMR may 
depend on the detailed band structure of the system, its 
angular dependence is essentially determined by the sym- 
metry properties of the SOC field. Here we investigate 
how the specific form of the TAMR angular dependence 
emerges from the properties of the SOC field. 

A phenomenological model which incorporate the ef- 
fects of the interference of Bychkov-Rashba and Dressel- 
haus SOCs was recently developed to explain the in-planc 
TAMR in (001) ferromagnet /semiconductor/normal- 
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metal (F/S/NM) MTJs,^iiiii^ In particular, it was shown 
that, in spite of its relative simplicity, the model was able 
to reproduce the two-fold symmetric angular dependence 
of the in-plane TAMR experimentally observed in (001) 
Fe/GaAs/Au MTJsj^i^i In such heterojunctions aU the 
involved materials are cubic in their bulk form. There- 
fore, the two-fold anisotropy of the in-plane TAMR must 
originate from the interfaces. Here we generalize the 
model and provide a unified qualitative description of 
the angular dependence of both the in-plane and out- 
of-plane TAMRs in (001), (110), and (111) MTJs. We 
consider systems in which the SOC originates from struc- 
ture inversion asymmetry (Bychkov-Rashba-like SOC) 
and/or bulk inversion asymmetry (Dresselhaus-like SOC) 
and predict new forms of the TAMR angular dependence 
which could be tested in future experiments. The effects 
of uniaxial strain arc also discussed. 



II. THEORETICAL MODEL 

We consider a MTJ composed of a ferromagnetic elec- 
trode and a normal metal counter-electrode separated by 
an insulator or a semiconductor barrier. However, our 
conclusions are also valid for the case of MTJs with two 
ferromagnetic electrodes whose magnetizations are par- 
allel to each other, since such systems are qualitatively 
similar to the case of MTJs with a single ferromagnetic 
lead. 

The z direction is fixed along the normal to the fer- 
romagnetic layer (i.e., parallel to the growth direction). 
The effective spin-orbit interaction corresponding to the 
nth band can be written as 

i?so = w„(k) ■ cr, (3) 

where w„(k) = [i(;„a;(k), i(;„j^ (k), w„2(k)] is the effective 
SOC field associated to the nth, k is the wave vector, 
and cr is a vector whose components arc the Pauli ma- 
trices. Equation ([3]) is quite general, since by now we 
have not considered any specific form for the SOC field. 
The detailed form of the SOC field can be quite compli- 
cated as one goes away from the center of the Brillouin 
zone and quite different from band to band, as recently 
demonstrated by first-principles calculations^^. 

Due to the presence of the spin-orbit interaction, the 
transmissivity T„(k, m) corresponding to the nth band 
becomes dependent on the magnetization direction m. 
Assuming that the strength of the SOC field is small rela- 
tive to both the Fermi energy and the exchange splitting, 
one can expand the transmissivity in powers of w„(k). 
For a given n and k there are only two preferential di- 
rections in the system, defined by rh and w„. Since the 
transmissivity is a scalar function, it can be written, to 
second order in the SOC field strength, in the for m^i^^i^^ 

r„(k, m) « al°j(k) + al'j(k)[m • w„(k)] + 
al?,)(k)|w„(k)p + 4^J(k)[m . w„(k)]^ (4) 



which represents the most general expansion (up to sec- 
ond order) of a scalar function (the transmissivity) in 
terms of two vectors (m and w„). Note that the argu- 
ments used for obtaining Eq. ^ are also valid for MTJs 
with two ferromagnetic electrodes whose magnetizations 
are parallel to each other along the direction m. The 
expansion coefficients cl-^^ (i = 1,2; j = 0, 1,2) refer to 
the system in the absence of the SOC field and therefore 
do not depend on lii. Since these coefficients reflect the 
cubic symmetry of the involved bulk materials they obey 
the relations a'fj{k^,ky) = a^^^-k^, -fc^), af^(.^x, ky) = 
o-tni-^x, ky), and a[l^{k^,ky) = a'fj^ky.k^). Cases in 
which the involved materials have other than cubic sym- 
metry in their bulk form can be treated analogously. 

Within linear response theory, the conductance G 
through the MTJ is determined by the states at the Fermi 
energy Ep and the k dependence reduces to the in-plane 
k|| dependence at E ~ Ep- One can then write 

G('ii)-|^E(^"(k|hrii)), (5) 

n 

where = 2e^ /h is the conductance quantum and (...) 
denotes evaluation aX E = E^ and integration over k|| . 

The time reversal symmetry implies that T„ (k, rh) = 
r„(— k, — ih) and w„(k) = — w„(— k). It follows then 
from Eq. ([¥]) that the first-order term in the expansion 
must be an odd function of k and will, therefore, vanish 
after integration over k|| . As a result the conductance 
can be rewritten as 

G(m) = GW+GL^<]+Gil„(m), (6) 

where G^°-' is the conductance in the absence of SOC, 
GL^,'cx(a(^,'(k|||w„(k||)|2)) and 

Gilo(m) - |^^(4'i(k||)[m. w„(k||)]2), (7) 

n 

are the isotropic and anisotropic SOC contributions, re- 
spectively. In terms of the components of m and the 
SOC field, Eq. reduces to 

Ginlo(^,'^) = |^Tr[AM(0,0)], (8) 

where A and M(0, (p) are matrices whose elements are 
given by 

^iO = X]^"2n(k||)WmW,y) (i, j = x,y, z) (9) 
n 

and 

My(0,<^) = m,(0,0)mj(0,0) (i, j = x,y,z), (10) 
respectively. 

Equations ([5|)- ([nH) are quite general and reveal how the 
symmetry of the SOC field can lead to the anisotropy of 
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the conductance. Further simphfications of these expres- 
sions can be realized by taking into account the proper- 
ties of the specific form of the SOC field. In what follows 
we focus on the case in which the components of the SOC 
field can be written as Wni = d„i • k|| . As shown bellow, 
many relevant physical situations correspond to such a 
case. The matrix elements in Eq. then reduce to 

Ay = ^ C„ {dni ■ dnj) . (11) 

n 

In obtaining Eq. pip we took into account the four- 
fold symmetry of the expansion coefficients (k|| ) from 
which follows that the only non- vanishing averages are of 
the form c„ (a^^j (kjO'fc^) = (4^^ (k,, )fc2). By using the 
Eqs. ([H), Ho]), and (fTl|) the anisotropic part of the con- 
ductance can be rewritten as 

Ginlo = 1^ E <^)™.(^: ^) (dr»: • d„,) . (12) 

The dependence of the TAMR ratio on the magneti- 
zation direction is determined by the anisotropic part of 
the conductance. Thus, Eq. (|12p is our starting formula 
for discussing important particular cases. 

III. RESULTS 

We first neglect the effects of strain and focus on the 
particularly relevant case in which the SOC field results 
from the interference of the Bychkov-Rashba and Drcs- 
sclhaus SOCs. Later on we shall consider also MTJs with 
SOC induced by uniaxial strain. 

The Bychkov-Rashba SOC originates from the struc- 
ture inversion asymmetry (SIA) of the junction and is 
basically determined by the strong electric fields at the 
interfaces of the tunneling barrier. It is present, for ex- 
ample, in MTJs with the left and right electrodes made 
of different materials and therefore with broken inver- 
sion symmetry. The Dresselhaus SOC results from the 
bulk inversion asymmetry (BIA) of one or more of the 
constituent materials. Typical materials with BIA are 
the zinc blende semiconductors. Thus, the Dresselhaus 
SOC can be relevant for MTJs with non-centrosymmetric 
semiconductor barriers. 

The specific form of the SOC field depends on the 
growth direction of the heterostructure. Below we an- 
alyze the most relevant cases, corresponding to MTJs 
grown in the [001], [110], and [111] crystallographic di- 
rections. 

A. (001) MTJs with axes x || [110], y || [110], z || [001] 

In this case the SOC corresponding to the nth band 
containing both Bychkov-Rashba and Dresselhaus terms 
is given byiii^i^ 

HsO = ("n - ln)KCy - (a„ + Jn)kyCrx, (13) 



where a„ and 7„ are the corresponding Bychkov-Rashba 
and Dresselhaus parameters, respectively. 

One can extract the components of the SOC field by 
comparing Eqs.© and (|13p . It follows then from equa- 
tion (fT2)) that the angular dependence of the anisotropic 
conductance is given by 

Gilo = E ^" + ll) + 2«„7« cos(20)] . 

n 

(14) 

The expression above together with Eqs. H]) and ^ lead 
to the relations corresponding to case A in Table 
The obtained TAMR coefficients, which are valid up to 
second order in the SOC field, reveal a clear distinction 
between the in-plane and out-of-plane configurations in 
[001] MTJs: while for a finite out-of-plane TAMR the 
presence of only one of the SOCs suffices (i.e., it is suf- 
ficient to have ^ or 7„ 7^ 0), the two- fold symmet- 
ric in-plane TAMR appears because of the interference 
of non- vanishing Bychkov-Rashba and Dresselhaus SOCs 
(i.e., both a„ and 7„ have to be finite).— This explains 
why a finite out-of-plane TAMR appears in MTJs such 
as Fe(001)/vacuum/Cu(001) in which only the Bychkov- 
Rashba SOC is present^ and is in agreement with the re- 
cent observation of the in-plane TAMR in epitaxial (001) 
Fe/GaAs/Au MTJs?^ where due to the presence of the 
non-centrosymmetric zinc blende semiconductor GaAs as 
the barrier material not only the Bychkov-Rashba but 
also the Dresselhaus SOC become relevant. In both the 
in-plane and out-of-plane configurations, angular depen- 
dencies of the form TAMRj^'io] {<P) oc [1 - cos(2(^)] and 
TAMRfi"!*;,] (9) cx [cos(26') - 1] (see Table H]) have been 
experimentally measured 41^ 1^°'^^ 

The results displayed in Table |T] suggest the possibil- 
ity of using different configurations and reference axes as 
complementary setups for TAMR measurements. In par- 
ticular, our theoretical model predicts that in the regime 
an ~ 7n the out-of-plane TAMR with reference axis in 
the [110] is suppressed, while it remains finite if the com- 
plementary axis [110] is used as a referencei^ The oppo- 
site behavior, i.e., TAMR|\%] ^ and TAMRfj^i'o] = 0, 
is expected when a„ « — 7n. Another relevant regime 
occurs when a„ w for which the in-plane TAMR 
is expected to vanish (see case A in Table The 
existence of such a regime was previously invoked in 
Refs. [Mill for explaining the suppression of the in-plane 
TAMR experimentally observed in (001) Fe/GaAs/Au 
MTJsj^ Our theory predicts that although the in-plane 
TAMR vanishes, the out-of-plane TAMR should remain 
finite in such a regime. In fact, in the regime Q!„ « 
the amplitude of the out-of-plane TAMR constitutes a 
direct measurement of the effects of BIA in the non- 
centrosymmetric barrier. 

By combining the results shown in Table |T] one can find 
expressions such as 

TAMRi^io](90°) = TAMR|\"i*o](90°) ~ TAMRfi'\*o](90°), 

(15) 
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TABLE I: TAMR coefficients in units of .g(,/[167r^(G'°' + G\^1)] for different structures, reference axes [x], and configurations. 



Case 


Structure 


[x] 


in-plane TAMRj^', ((?i) 


out-of-plane TAMR^f (S) 


A 


(001) MTJ 


[110] 


4[1 - cos(2<;/))] J2„ c„a„7„ 


[cos(26i)-l]E„c„(a„+7n)^ 






[110] 


4[cos(2(/)) - 1] J2„ Cnan'yn 


[cos(2e)-l]E„c„(a„-7„)2 


ri 


(iiOj Mi J 


[110] 


[1 - cos(2(^)J 2^,^ c„(q„ - /i„) 


[cos(2f^) - 1] c„(q!„ - A„) 






[001] 


[cos(20) - 1]E„c„(q^ _/32) 


2 E„ c„ [A^ - (/3„ sin e + X„cosef] 


c 


(111) MTJ 


[112] 





[cos(26i)-l] ^^c„(a„+7n)' 






[110] 





[cos(26l)-l]E„c„(a„+7„)" 


D 


strained (001) MTJ 


[110] 


4[1 - cos(2(^)] c„rj„fj.„ 


[COS(20) - 1] J2nC"{Vn+ fJ-nf 






[110] 


4[cos(2(^) - 1] c„ri,^firi 


[cos(2e)-l] E„c„(??„-/in)' 



SOC resulting from the BIA. Note that because of the 
reduced symmetry of the (110) structures with respect 
to the (001) MTJs, in the present case the usual SIA- 
induced SOC acquires, in addition to the usual Bychkov- 
Rashba SOC, an extra contribution which leads to q:„ ^ 
fin in Eq. 



Proceeding in the same way as in Sec. IIII Al we obtain 
the following relation for the anisotropic contribution to 
the conductance 



which correlate the in-plane and out-of-plane TAMR co- 
efficients and can be experimentally tested. 

B. (110) MTJs with axes x || [110], y |1 [001], z [| [110] 
In this case the SOC is described hy^^i^ 

Hso = ankyCJx - PnkxO-y - Xnkx<7z- (16) 

Here a„ and /?„ are parameters related to the SIA- 
induced SOC, while A„ characterizes the strength of the 



G 



sin^ 6' + cos^ 6* -f /3„ A„ sin sin(26')] 



(17) 



The corresponding TAMR coefficients are given in Table 
m The relations for case B show that the angular depen- 
dence of the TAMR in both the in-plane and out-of-plane 
are similar to the ones obtained for the (001) MTJs [com- 
pare the cases A and B in Table|l]. However, their physi- 
cal origin is now different. In the present case the in-plane 
TAMR originates from the SIA induced SOC while the 
out-of-plane TAMR has contributions arising from both 
SIA- and BIA-like SOCs. Thus, our model predicts that 
in (110) MTJs it could be possible to observe the TAMR 
in the two configurations even if the tunneling barrier is 
composed of a centrosymmetric material. Another ob- 
servation is that the out-of-plane TAMR with reference 
axis along the [IlO] direction could be suppressed if under 
some given conditions the regime q;„ = ±A„ (for the rele- 
vant to transport bands) is realized [see the out-of-plane 
TAMR in case B of Table |l]. In such a case, however, 
the out-of-plane TAMR with [001] as the reference axis 
should remain finite. 



C. (Ill) MTJs with axes x |[ [112], y [| [110], z [] [111] 
For this structure the SOC is given by^^ 

^^SO = {O-n + ln)ikx(Ty ~ kyCJx), (18) 

where a„ and 7„ are the parameters characterizing the 
strengths of the SIA- and BIA-like SOCs, respectively. 

After computing the anisotropic part of the conduc- 
tance we obtain 

Gilo = ^E^"K+7„f. (19) 

n 

This relation leads to the TAMR coefficients given in 
Table HI for the case C). 

In the present case the prediction of a vanishing in- 
plane TAMR is remarkable. We have checked that even 
if the cubic in k terms are included in the SOC field, 
the in-plane TAMR still vanishes. This could be used 
for experimentally exploring the origin of the TAMR. If 
a suppression of the in-plane TAMR is experimentally 
observed in (111) MTJs, it will be a strong indication 
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that indeed the mechanism behind the TAMR is the SIA- 
and/or the BIA-hke SOCs. On the contrary, if no sup- 
pression of the in-plane TAMR is observed, the role of 
these spin-orbit interactions as the origin of the TAMR 
can be questioned. 

Another interesting issue is the possibility of reaching 
the condition a„ = — 7„ (for the bands relevant to trans- 
port), which leads to a vanishing out-of-plane TAMR 
(if only the linear in k terms in the SOC field are rel- 
evant), in addition to the above discussed suppression of 
the TAMR in the in-plane configuration. 



where Uij are the components of the strain tensor and 
and 7„ are materials parameters. The SOC in Eq. (j20p is 
quite rich and suggests the possibility of engineering the 
strain (see, for example, Ref. [28l ) in order to manipulate 
the behavior of the TAMR. Here we do not consider all 
the possibilities but focus, for the sake of illustration, on 
the case of an in-plane uniaxial strain such that the only 
non- vanishing components of the strain tensor are u^x ~ 
Uyy 7^ Uxy = Uyx- Thc cxistcnce of a similar strain was 
initially assumed for explaining TAMR experiments in 
(Ga,Mn)As/A10x/Au MTJs-ii^ For the in-plane uniaxial 
strain Eq. (^0]) reduces to 

^SO = rjnikxCTy - kyCFx) + fJ-nikxC^x - kyCTy), (21) 

where we have introduced the strain-renormalized SIA 
and BIA parameters ?7„ = Q:„Uxy and /i„ = "JnUxx, re- 
spectively. The corresponding anisotropic contribution 
to the conductance is then given by 

^inlo = ^""^^ E M + + 2^"M« sin(2<^)] . 

n 

(22) 

We note that Eq. ((2T|) has the form of interfering 
Bychkov-Rashba and Dresselhaus SOCs in (001) struc- 
tures with X II [100] Jii^ Therefore, Eqs. ^ and ^ 
are similar. The angle in Eq. (|14[) is measured with 
respect to the crystallographic direction [110] while in 
Eq. (12211 it is defined with respect to the [100] axis. Thus, 
by making the transformation — > (/> + 7r/4 in Eq. (j22p 
one recovers a relation similar to Eq. (jl4p . Consequently, 
assuming the direction [110] as the reference axis for mea- 
suring the magnetization direction the results for the in- 
plane and out-of-plane TAMR coefficients in (001) MTJs 
with in-plane uniaxial strain [see case D in Table |T] are 
essentially the same as in the case discussed in llll Al but 
with renormalized spin-orbit parameters which now ac- 
count for thc strain effects. 



D. Uniaxial strain in (001) MTJs with axes 

X II [100], y II [010], z II [001] 



In our previous analysis we have disregarded the ef- 
fects of the strain induced SOC which could be relevant 
for structures whose constituent materials have a sizable 
mismatch in their lattice constants. For a (001) MTJ the 
SOC induced by strain is, in general, given by22 



(20) 

I 

In our investigation we have assumed specific well- 
known forms for the SOC field. For some systems, 
however, the form of the SOC field may not be a pri- 
ori known. In such a case one could use Eqs. l|5))- pUl) 
(which are general) and contrast them with complemen- 
tary TAMR measurements in both the in-plane and out- 
of-plane configurations in order to deduce the symmetry 
properties of the SOC fields. 

All the calculated TAMR coefficients, if not zero, 
show a two-fold symmetry in the {6, (^)-space, which 
is the symmetry that has been observed in the 
experiments4i^'^ '^°i^^ Our results are valid up to the sec- 
ond order in the SOC field. In particular, our predic- 
tions for vanishing TAMR under certain conditions, may 
change when higher orders in the SOC field become rele- 
vant. The next higher order contributions in the expan- 
sion in Eq. ^ which do not vanish after averaging are 
those containing the fourth order in the SOC field terms 
and terms of the fourth order in the cosine directions of 
m which describe the four-fold symmetry inherent to the 
involved bulk ferromagnet. These fourth order terms lead 
to four-fold symmetric corrections to the TAMR which 
may be finite even for (001) MTJs with centrosymmetric 
barriers for which the second order in-plane TAMR cal- 
culated here vanishes. Although the two-fold character of 
the TAMR is, in general unchanged by these corrections, 
they may influence the shape of its angular dependence. 
Thus, for the kind of systems here considered (see Ta- 
ble [l| any deviation from the 8 — like polar shape of 
the TAMR [sec, for example. Fig. 2 in Ref. Q is inter- 
preted in our theory as a manifestation of higher order 
contributions and/or strain effects Deviations from the 
8 — like polar shape of the TAMR have been experimen- 
tally observedi^iiS In fact, it has been shown that these 
deviations may appear by increasing the bias voltage;^ 
which within the present approach can be seen as an in- 
dication of higher order in the SOC field terms turning 
relevant at sufficiently high bias. 



I 

Hso = an [{uzxkz - Uxyky)ax + [uxykx ~ Uyzkz)ay + [uy^ky - Uzxkx)<Jz] 

+ Tn [kx{Uyy 'Uzz}'^x "I" ky(Uzz Uxx)o'y + kziUxx Uyy)a z\ j 
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In all the above discussions small magnetic fields with 
negligible orbital effects were assumed. It has recently 
been observed in Fe/GaAs/Au MTJs that for high mag- 
netic fields the orbital effects do influence the in-plane 
TAMR.^ A version of the phenological model here pre- 
sented which incorporates the orbital effects has recently 
been developed to qualitatively explain the magnetic field 
dependence of the in-plane TAMR experimentally ob- 
served in (001) Fe/GaAs/Au MTJs,^ 

IV. SUMMARY 

We formulated a theoretical model in which the way 
the TAMR depends on the magnetization orientation of 



the ferromagnetic electrode in MTJs is determined by the 
specific form and symmetry properties of the interface- 
induced SOC field. By using the proposed model, we 
deduced the angular dependence of the TAMR for various 
systems in dependence of their symmetry under spatial 
inversion and their growth direction. The effects of in- 
plane uniaxial strain were also investigated. 
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